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Abstract. We consider nonlinear Sclirodinger equations witli eitlier local or nonlocal 
nonlinearities. In addition, we include periodic potentials as used, for example, in matter 
wave experiments in optical lattices. By considering the corresponding semiclassical scal- 
ing regime, we construct asymptotic solutions, which are concentrated both in space and in 
frequency around the effective semiclassical phase-space flow induced by Bloch's spectral 
problem. The dynamics of these generalized coherent states is governed by a nonlinear 
Schrodinger model with effective mass. In the case of nonlocal nonlinearities we establish 
a novel averaging type result in the critical case. 



1. Introduction 

Coherent states have been originally introduced in quantum mechanics to describe wave 
packets minimizing the uncertainty principle. This property makes coherent states highly 
attractive for the study of semiclassical asymptotics, see, e.g., lfT7llT5l[2TI . Indeed, it can be 
shown that for Schrodinger equations with sub-quadratic potentials, coherent states retain 
their shape, providing minimum uncertainty at all time in the quadratic case ifTSl . and up 
to Ehrenfest time in general 16 |. Recently, extensions to weakly nonlinear situations have 
been studied in |8 , 7 |. In addition, the semiclassical dynamics of coherent states under the 
influence of (highly oscillatory) periodic potentials has been investigated by the authors 
in IfTOl . In the present work we combine the effects coming from periodic and nonlinear 
potentials. 

To this end, we consider nonlinear Schrodinger equations which, after scaUng into 
dimensionless coordinates, appear in the following semiclassical form: 

2 

where t e M, a; e W\ and o? e N denotes the spatial dimension (usually d = 3). Moreover, 
£ e (0, 1] denotes a (small) semiclassical parameter, i.e. a dimensionless rescaled Planck's 
constant. The factor e" measures the (asymptotic) strength of the nonlinearity: the larger 
the a > 0, the weaker the nonlinear effects. In the following, we shall allow for two 
different types of gauge invariant nonlinearities: 

• Local nonlinearities: /(l^/'^P) — A|7/'^p'^, with cr G N and A G M, allowing for 
focusing (attractive) and defocusing (repulsive) situations. 

• Nonlocal nonlinearities of convolution type: — ji/i^p, with i^(a;) G 
R a given interaction kernel. 

Finally, the term V-pcaix/e) denotes a highly oscillatory periodic potential. More precisely, 
let r ~ Z'' be some regular lattice, then we assume that for all y G M"^': Vpcr(y + 7) = 
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Vpci-{y) with 7 € r. In addition, we shall assume Vpcr G C°°(M'^). Equation jl.lb 
describes the propagation of waves on macroscopic length- and time-scales, i.e. over many 
periods of the periodic potential. The parameter e ^ 1 consequently describes the ratio 
between microscopic (quantum mechanical) and the macroscopic scales. 

Nonlinear Schrodinger equations with periodic potentials arise in various physical con- 
texts: A by now classical example is the mean-field description of electrons propagat- 
ing in a crystalline solid [ Id, 26J under the additional influence of a self-consistent elec- 
tric field. The latter is usually modeled by means of a nonlocal Hartree nonlinearity 
= I'^'^P * 1/1 ■ I' see, e.g., |T| for a semiclassical study via Wigner measures. 
Another situation in which dl.lb applies is the description of Bose-Einstein condensates in 
so-called optical lattices, cf. lfT2l . In the regime of dilute gases, such condensates can be 
modeled by the Gross-Pitaevskii equation with cubic nonlinearity (7=1, cf. ^TI\. Note 
however, that other nonlinearities also arise, see, for example, [5| where a nonlocal term 
is used for the description of superfluid Helium. In addition, strong magnetic confinement 
allows for the experimental realization of quasi-one-dimensional (cigar-shaped) conden- 
sates, or quasi two-dimensional condensates, motivating the fact that we consider ( II. lb in 
general dimensions d e N, cf . ifTSll . A third example for the appearances of ( 11.11 ) stems 
from the description of wave packets propagating within nonlinear photonic crystals |24| 
and where the nonlinear response of the media is modeled via a Kerr nonlinearity a = 1. 
In this case, the underlying assumption in the derivation of (11.11 ) is the existence of a pre- 
ferred direction of propagation, implying that the appropriate model is stated in dimension 
d= 1. 

In all of these situations, the joint effects of nonlinearity, periodicity and dispersion (or, 
quantum pressure), can lead to the existence of stable localized states conserving the form 
upon propagation and collisions. Gap solitons, discrete breathers and compactons are ex- 
amples of such states. Here we shall present another possibility, which will arise from 
semiclassical description via coherent states. Even though the nonlinearity in our case is 
weaker than in the above mentioned situations (due to the fact that a > 0), the obtained as- 
ymptotic solutions will nonetheless experience nonlinear effects in leading order, provided 
a is of critical size (to be made precise later on). The latter will depend on the precise form 
of the nonlinearity. 

To present our results, we recall the classical Block eigenvalue problem ll28l : 

(1.2) H(k)xm{;k)= E„,{k)xm{-,k), m e N. 

Denoting by Y the centered fundamental domain of F, Em{k) G R and Xm(-, k) denote, 
respectively, the rn-th eigenvalue/eigenvector pair of 

(1.3) 77(fc) = i (-zV, + fc)' + 1/pcr (y) , yeY, 
parametrized by the crystal momentum fc e ~ T''. We shall assume that at i = 0, 

(1.4) ^-^e-^"'u, (^) ,Po) e^P-(--'«)/% 

where mq denotes some smooth and rapidly decaying profile. In other words, the initial 
data -00 can be approximated by a highly oscillatory Bloch eigenfunction Xm modulated 
by a (generalized) coherent state, i.e. a wave function which is localized both in space and 
in frequency. 

Remark 1.1. In particular, the choice uq{z) = exp(— |2:p/2) yields a classical coherent 
state, i.e. ground state of the harmonic oscillator potential, as modulation. The states ( 11.4b 
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are more general, though, since we can allow for any uq E S{M.'^), the Schwartz space of 
rapidly decaying, smooth functions. We also remark that the same class of initial data have 
recently been considered in |TOl, where the situation of linear Schrodinger equations with 
combined periodic and slowly varying external potentials has been considered (see iB] for 
more details). 

Provided such initial data, we shall show that the solution of ( II. Il l can be approximately 
(in a sense to be made precise) described by the following semiclassical wave packet: 

(1.5) r{t,x)^^^^^{t,x) ■.^e-'^l'u (t,^^^) X™ (f ,Po) e^^"(*'-)/- 

where q{t) = qo + f^kEmiPo) describes the macroscopic shift of the centre of mass and 
the highly oscillatory phase function is 

(1.6) 4>m{t,x) =po ■ {x- go) - tE,n{po)- 

To this end, we need to give sense to the group velocity VkE„i{k) and thus, we have to 
assume from now on, that: 

(1.7) E,n{k) is a simple eigenvalue in the vicinity of k ~ pq. 

(Of course, po G has to be understood modulo T* in this case.) In other words, we 
have to avoid that two Bloch bands cross at po, that is Emipo) = En{po), for m ^ n. It 
is known that at such crossing points E,n{k) is no longer differentiable, causing the above 
asymptotic description (which is based on an adiabatic decoupling of the slow and fast 
degrees of freedom) to break down. 

Remark 1 .2. Clearly, the non-crossing condition given above restricts our choice for the 
initial wave vectors po £ Mf^. It is known however that the set of band crossings has 
Lebesgue measure zero. For example, in the case d = 1, band crossings can only occur at 
fc = or at the boundary of the Brillouin zone. 

So far, we have not said what determines the profile u — u{t, z) appearing in (ll.5l l. Its 
time-evolution depends on the strength of the nonlinearity, i.e. on the size of a > in the 
case of local nonlinearities (the situation for nonlocal ones will be described later on). We 
shall find that the critical size is ac = I + da/2 and when a = ac, u solves the following 
homogenized nonlinear Schrodinger equation 

(1.8) + ^div^ ((V^-Em (po)) • V^) M = A„|u|^'^u, U|t=o = ""o, 
with effective coupling constant 

Am = \ Ixm {y,po)f"^'^ dy. 

Note that the dispersive properties of ( 11.81) are determined by an effective mass matrix 
Wf.Em (po) G M.'^^'^, which itself depends on the choice of the initial wave vector (cf. lfT4l 
for a recent study). 

In the next section, we derive this effective mass equation from multi-scale expansion. 
A rigorous stability result is then proved in Section [3] The case of nonlocal nonlinearities 
is treated in Section ID As we shall see, in situations where the kernel K is homogeneous, 
the critical value ac depends on the degree of homogeneity (like in the case of a local non- 
linearity), and the analogue of ( II.8I 1 is an envelope equation with a nonlocal nonlinearity. 
On the other hand, if the kernel X is a smooth function, then etc = 1, and in sharp contrast 
with the other situations studied in this paper, the analogue of ( II.8I 1 for a ~ ac is, found to 
be a linear equation, see Section |4~2| 
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2. Multi-scale expansion 

2.1. The hierarchy of equations. Except for the treatment of the nonlinear term, we re- 
sume the strategy followed in ifTOl . We seek the solution ip'^ of (II. lb in the form 



(2.1) 



i,%t, x) = e-^/^W {t, ^^1^, f ) e"^'"(*■-)/^ 



where the phase (/)m(t, x) is given by (|1.6l l and the function W ~ l^^{t, z, y) admits an 
asymptotic expansion 

(2.2) 

with smooth profiles Uj which, in addition, are assumed to be F-periodic with respect to 
y. For ifj'^ given by the ansatz ( 12.1b . we compute 



ledtr + Y - ^pcr (-) = e-^'^e^*-'^- ^ e='^{t, z, y) 
with 

bl = ~iq{t) ■ VM' + (V,y • + ipQ ■ V,W, 

bl = idtW + ^AM'. 

Using ( |1.6b and the fact that q{t) = W kEmipo), we can rewrite 6q, &f as 

bl =l{pQ~ VkErn (PO)) ■ ^zU' + (Vy • V,)W^ 

where H{po) is the Bloch Hamiltonian ( 11.3b evaluated at A: = po (again, this has to be 
understood modulo F*). Introducing the following linear operators: 



Lq = E„i {pa) - H (po) , Li=i{po - VkE„i (po)) • + • Vz, L2 = i9t + -A^ 



and expanding W in powers of e G (0, 1], we consequently need to solve the following 
hierarchy of equations: 

r LoUo = 0, 
(2.3) J LoUi+LiUo = 0, 

[ L0U2 + LiUi + L2U0 = F{Uo), 

where, for 0;^ = 1 + d(T/2, we find: 

{0 if a > ttc, 
MUo^Uo if a = 

In the next subsection, we shall focus on the resolution of ( 12.31 ). 
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2.2. The effective mass equation. Given the form of Lq, the equation LqUq — implies 

(2.4) UQ(t,z,y) = u{t,z)xm{y,Po) ■ 

By Fredholm's alternative, a necessary and sufficient condition to solve the equation LoUi + 
LiUq = 0, is that LiC/q is orthogonal to kcrio^ that is: 

(2.5) (Xm, LlUo) i2(y) = 0. 

Given the expression of Li and the formula (I2.4l l. we compute 

LiUo = i{po- VkErn (Pq)) ■ ^zU{t, z)Xm (V , Po) + '^yXm iv , Po) ' V^u(i, z). 

Now, we make use of the algebraic identities derived in Section |A] In view of ( IA.2l i. we 
infer that (12.51) is automatically fulfilled. We thus obtain 

Ui{t,z,y) = ui{t,z)x,n {y,po) + ui{t,z,y), 

where uj^, the part of Ui which is orthogonal to ker Lq, is obtained by inverting an elliptic 
operator: uf- = —Lq^LiUo- Note that the formula for LiUq can also be written as 

LiUo = -iVfc {Ejn (pq) - H{po))xm {y,Po) ■ V^uit.z). 

Taking into account ( lA.ll l. this yields: u^{t, z, y) — ~iV kXm {y,Po) • '^z'u{t, z). At this 
stage, we shall, for simplicity, choose ui = 0, in which case Ui becomes simply a function 
of u: 

(2.6) Ui{t,z,y) = -iVkXm (2/,Po) • '^zu{t,z). 

As a next step in the formal analysis, we must solve L0U2 + LiUi + L2UQ = F{Ua). By 
the same argument as before, we require 

(2.7) (x™, LiUi + L2U0 - F{Uo))l2(^y) = 0- 
With the expression (12.61 1. we compute 

d 

LiUi = ^ (^{po - VkErn {pQ))j dktXm {v,Pq) " idy-k^Xm (y,Po)) d^.^^u, 



and we also have 



L2U0 ^ [idt + ^A^ ) u{t, z) Xm (y,Po) 



We consequently infer 

1 
2 



(Xm, LiUi + L2U0) ^2(Y) = idtu + -A^u 



Xrn) ^ ^ZjZi^- 

In the case a = etc., by making the last sum symmetric with respect to j and I, and using 
( IA.3b . we finally obtain the nonlinear Schrodinger equation dl.SI l with effective mass tensor 
M = V|£'m(po) and coupling constant 

In addition, we can write 

(2.8) U2{t,z,y) = U2{t,z)x,n {y,po) + U2{t,z,y), 

where = —Lq^ {LiUi + L2U0) . We shall also impose U2 = and thus U2 = U2. 
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Assume for the moment that we can solve ( 11.81) . Then, we have the following result, 
which establishes some basic regularity properties of our multi-scale expansion (where we 
denote by H'' the usual L^(R'^) based Sobolev space). 

Lemma 2.1. Suppose ( 11.71 ) holds true and let u G C([0, T]; H^), be a solution of ( 11.81 ) up 

to some T > 0. Then Uj £ C([0, T];H^-^ x W'^'°°{Y)),for j = 0, 1, 2. 

Proof. First note that {y,k) i— Xm{y,k) is smooth and bounded together with all its 
derivatives, provided il.li holds true. Having this in mind, the proof follows directly from 
the construction of {[/j}j=o,i,2 as solutions to the system ( 12.3b . □ 

Remark 2.2. Note that in the case a > ac nonlinear effects are absent at leading order 
since we obtain, instead of ( II.8I 1, a linear effective mass equation: 



This type of equation has been derived in |I2]|25|, using a different asymptotic scaling. 



In this section we shall make the computations given above rigorous and prove a nonlin- 
ear stability result. As a first step we need to guarantee the existence of a smooth solution 
to ( II.8I 1, at least locally in-time. 

3.1. Construction of an approximate solution. The dispersion relation of ( |1.8t is given 
by a real-valued symmetric matrix. Standard techniques (see, e.g., fTT\) yield the existence 
of a unique local solution, provided that the initial datum is sufficiently smooth: 

Lemma 3.1. Let uq E with k > d/2. There exists Tc G (0, +00] and a unique 
maximal solution u G C{[Q^Tc)\ H*') to ( 11.81 ), such that \\u{t,-)\\i^2 — ||uo||/^2. The 
solution is maximal in the sense that ifT^ < 00, then 



The solution u{t, •) may not exist for all times, even if A ^ 0, i.e. even if the nonlinearity 
in the original equation ( II. lb is defocusing. However, we can claim Tc — 00 in either of 
the following cases (see e.g. IfTTII ): 

• '^l.Emipo) is positive definite and X„i ^ 0, or 

• V^ii^m(po) is negative definite and Am ^ 0. 

On the contrary, if for instance V|£'m(po) is positive definite and Am < (focusing 
nonlinearity), finite time blow up (that is, Tc < 00) may occur, see, e.g., IIZTHTTI . This is 
the case typically if the initial datum is "too large": for any fixed profile uq £ iS(]R''), if 
one considers U|t=o = Auq, there exists Aq > such that for all A ^ Aq, Tc < oo. Note 
that in other situations, where the signature of V^i?m(po) is non-trivial (hence d ^ 2), the 
issue of global existence vs. finite time blow-up is an open question. 

Remark 3.2. Clearly, for a > ac these issues do not occur, since the leading order profile 
Miin solves the linear equation ( 12. 9b and hence exists for all times, Tc = +00. 

Lemma [TTI provides the existence of a local-in-time solution u of the effective mass 
equation. In view of the multi-scale expansion given in Section |2l we can thus define an 
approximate solution by 




3. Main results 



lim \\u{t, •)||^fc = +00. 



(3.1) 



Vipp(i, x) := e-^/^ {Uo + V^Ui + eC/a) t, 




NONLINEAR COHERENT STATES IN PERIODIC POTENTIALS 



7 



which satisfies the original equation jl.lb up to some remainder terms: 

^2 



app I T^app 



^£3/2 {r\+rl){t,z,y) 



The remainder terms are given by 

r'i(t,z,y) = L2Ui{t,z,y),r2{t,z,y) = LiU2(t,z,y) + ^/eL2U2{t, z,y), 



and = e"|V'appP'^V'app if " > "c, while if a = etc = 1 + 



l+df7/2i |2<T„/,e 



app I V'app g(i/4^ 



Uolt, 



X — q(t) X 



Uolt, 



X — q(t) X 



This, together with the regularity result established in Lemma IZTI then directly yields the 
following proposition. 

Proposition 3.3. Assume (11.7b and let a ^ Uc = 1 + da/2. Then, we can find V'app such 
that: 

1. For all T G [0, Tc), V'lpp has a coherent state structure on [0, T], in the sense that there 
exists C independent of e G (0, 1] such that, with 93^ defined in ( 11.5b , 

fora = ac, sup \\%pl^^{t, ■) - ip" {t, ■)\\L2(^d) Cy/e, 

tG[0,T] 

fora>ac, sup ||'0app(^, ') - '/'fin(^: ■)llL2(Rd) < Cv^e, 

t6[0,T] 

where (/sfj^^ is f/ze approximate solution constructed from uun, solving ( |2.9b - 

2. The function solves dl.ll ) m/? fo a swa/Z error: 



^edti^lpp + yAV^^pp = Fper (^) V^pp + A£"|V: 



app I T-app 



■ ew 



where the remainder term satisfies: for all T > 0, with T < Tc in the case a = ac, 
there exists C > independent of e > such that 



„min(a — Qc ,1/2) 



sup Ww^it, •)llL2(R<i) C ■ 
tG[0,T] 



if a > ac, 
\/e if a — ac- 



Note that because of the factor e in front of the time derivative, it is natural to represent 
a small error term as e times as small term. 



3.2. Nonlinear stability. It remains to prove nonlinear stability of the approximate solu- 
tion constructed above. For the sake of simplicity we shall do so only for a — ac and 
d = 1. The (physically less interesting) case a > ac can be proved analogously and a pos- 
sible generalization to higher dimensions is indicated in Remark [331 below. For Eq > 0, 
set 

wrwni sup iwrh^ + wed^rh-), 

0<esj£o 

which is equivalent to the usual iJ^-norm for every (fixed) e > 0. The approach that we 
present is similar to the one followed in 19] : First, we need to construct a more accurate 
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approximate solution than the one stated in Proposition l3.3l Taking the asymptotic expan- 
sion presented in Section |2] one step further, we can gain a factor ^ in Proposition [33 
More precisely, we can construct ^/'app such that; 



(3.2) 
and 



sup 

tG[0,T] 



^eS.V'Ipp + ^AV^pp = ^per (^) V^pp + Ae^+^/^lV^prV^^pp + 



ew 



app ■ 2 

where the additional factor ^/e is reflected in the error estimate 

(3.3) sup \\w'{t,-)\\fj, ^ Ce. 

te[o.T] 

Note that in this case the corrector Ui is not the same for V'app> since unlike what we 



have done in ^ ^2.21 we can no longer assume ui = 0. Rather, ui now solves an evolution 
equation, which is essentially jl.SI l linearized about u, with a non-trivial source term (see 
^ for more details). Therefore, the estimate ( 13.2b must be expected to be sharp in general. 

Having constructed such an improved approximation V'app we can state the following 
nonlinear stability result: 

Theorem 3.4. Let d = 1, a = 1 + a/2, a ^ N, and Assumption (11.71) hold. In addition, 
suppose that the initial data satisfy: 



(3.4) 



% - V'app|t=0 



app|i— 



Oil) 



Let T £ [0, Tc). Then, there exists Eq = £o{T) such that for e £ (0, Eq], the solution of 
( II. Il l exists on [0, T]. Moreover, there exists C independent of e G (0, Eq] such that 

sup \\r{t,-)-V%t,-)\\L-iR)^CV^. 

te[o,T] 

where ip'^ is defined in ( ll.5l l. 

Proof. The scheme of the proof is the same as the proof of Theorem 4.5 in [9 1, so we shall 



only give the main steps. Fix T < Tc and let 77^ 
exact and the approximate solution. It satisfies 



V'^ — ipipp be the error between the 



,edtv' + i-Ar' = ^p„. (J) + Xe'+-^' (iV'l'V^ - IV'^ppPV^pp) 

with ||?7j^j^Q||L2 = 0{e), ||e9a;77j^j^g||L2 ~ 0{1) by assumption. From |23|, we have: 

Lemma 3.5 (Moser's lemma). Let R > 0, s £ N and F{z) = \z\'^" z, cr e N. Then there 
exists C = C{R, s, a) such that if v"^ satisfies 



and S"^ satisfies |l<5'^|l^o 



^ R, then 



^ i?, ^ /3 s, 
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We apply this lemma with v"^ = e^^'^ipl^p, and s = 0, s = 1 successively; there exists 
R> independent of e G (0, 1] such that 



sup V Ued^fv'it) 

*6[0,T] ^^o /l 



€ R. 



Set S'^ = E^l'^jf . By assumption and the Gagliardo-Nirenberg inequality, 

(3.5) W\,^^\\l^ = ei/^hf,.olU- ^ e^l'^2e-^l^\^,^Xl^\\ed^n\,^XL^ ^ Ce'/\ 

As long as ^ R, energy estimates and Moser's lemma with s = yield 



M{t)\\L2^Mm\L-+C U{s)\\L^ds+ W{s)\\L.ds, 

Jq Jq 
where we have used the homogeneity of F. By Gronwall's Lemma, for t ^T: 

Applying the operator edx to the equation satisfied by rj'^, we infer similarly 



^ \\w'{s)\\L2d^ sc Ce. 



^ Jo 



where the last term stems from the relation [edx, Vpci-{x/e)] = dyVpci{x / e) E L°°, since 
Fpcr is smooth and periodic. Thus, 

\\edxv'{t)\\L- ^C + C f \\edxV%s)\\L2ds + Ct. 



Gronwall's lemma now yields Wedxf]^ {t)\\L2 ^ C{T). In view of the Gagliardo-Nirenberg 
inequality. 

For e sufficiently small (depending of T), ||(5^(t)||L=° ^ i? for alH e [0, T], and the result 
follows from a bootstrap argument. □ 

The above theorem shows nonlinear stability of the approximate solution up to times 
of order 0(1), i.e. independent of e, provided that the initial data are well-prepared, in 
the sense given in ( |3.4| l. Essentially this means that i/)g contains not only C/q, but also Ui 
associated to V'app- We shall not insist further on this aspect, which is probably a technical 
artifact, and remark that in the linear case a stronger result is valid, see | lOJ where stability 
is proved up to the so-called Ehrenfest time 0(ln 1 /e), and no well-preparedness as in ( 13. 4t 
is needed (an initial error 0(e'') for some r > suffices). 

Remark 3.6. If a; G M'' with rf ^ 2, the proof can be easily adapted, provided an even 
better approximate solution is constructed. The reason is that, instead of ( 13.5b . one needs 
to rely on the following Gagliardo-Nirenberg inequality 

Ce-'^'^WW'-^ll^;^ II |£Vr<5^||f,f;i), for. > d/2. 

Thus, in order to account for the singular factor e^'^/^, one is forced to construct an ap- 
proximate solution V'app to a sufficiently high order in e (see |9| for more details). 
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4. The case of nonlocal nonlinearities 



In this section we shall show how to perform the same asymptotic analysis as before in 
the case of nonlocal nonlinearities. In other words, we consider 

2 

(4.1) iedtr + = ( r + e'^iK * \r?w 

with K{x) e M some given interaction kernel. In the following we shall focus on two 
particular choices of interaction kernels K which are physically relevant. 

4. 1 . Homogeneous kernels. In this subsection we shall consider functions of the form 

K{x) = A|a;|^, A G R, with ^ e M \ {0} such that - min(2, d) < fi^2. 

For example, the choice fi = —1 in d ^ 3 corresponds to the classical Hartree nonlinearity, 
modeling a self-consistent, repulsive (A > 0) Coulomb interaction. The case fi > has 
been recently studied in 1 19|. 

Like in the case of local nonlinearities, the critical exponent Qc depends on the homo- 
geneity /i, namely ac = 1 — /i/2. This can be seen as follows: We plug the ansatz 

X - g(t)~ 



Xr. 



into the convolution term (Nl^ * IV' 



U t, 



This yields 

2 

Xr 



-,Po 

e 



dC 



We want this term to be of order 0(e) in our asymptotic expansion, to mimic the approach 
presented in S|2] In this case, it will consequently appear within leading to the effec- 
tive mass equation. In order to show that this is indeed the case, we rewrite the initial 
convolution as 



A-fi/2-d/2 



U t 



Xr. 



dt 



and use the substitution z = {x — q{t)) / %/e in the envelope u, and y ~ xjem Xm- 



-l-/j/2-rf/2 



U \ t, Z — 



y - -,Po 

e 



Setting C = ^/y/s, this can be written as 



c 



Po 



dC 



\cnuit,z-o\' 

Then, the following averaging result can be proved: 

Proposition 4.1. Let (|1.7| i hold true and assume that ( t-^ \u {t, z ~ ()\^ is in 
Then, for all k € Y* , it holds 



\(:\'\u{t,z-o\' 



Xm y 



,k 



dC- 



\C\nu{t,z-0\ dC 



In addition, if C, i— )■ \u {t, z — C,)\ is in W'^'^(W^), then the above convergence holds 
with an error of order ©(^/e). 

This result can be seen as a variant of the two-scale convergence results introduced in 
||20l [n, and used in ||2l. The main difference here is the convolution structure. 



NONLINEAR COHERENT STATES IN PERIODIC POTENTIALS 



II 



Proof. We decompose y i-> |Xm(yjfc)P into its generalized Fourier series (recall that 
r ~ Z'') and write 

2 



icn^(^,^-c)r 



C 



7er 

By Riemann-Lebesgue lemma, for each term with 7^0, the limit, as e — 0, is zero. 
Then only the term corresponding to 7 = remains, with 



co(/c) 



IXm {v,k)\^ dy = 1, 



since the eigenfunctions Xm(-, k) form an orthonormal basis of L^(y). By the Dominated 
Convergence Theorem, we can exchange the sum over 7 € F and the limit e — > in the 
above computation provided that C, H> \u {t, z — C)l^ is in (with an error o(l)). 
In the case where the function is in W^'^ we obtain an error 0{y/e). The reason is that 
the coefficients (c-y)-ygr decrease rapidly for large I7I, since y i-> \xm{y, is smooth, 
provided Assumption ( 11.71 ) holds true and thus we can perform an integration by parts, and 
use dominated convergence again. □ 

Assuming that u is sufficiently smooth and decaying, we can use the above averaging 
result and perform the same asymptotic expansion as given in Section |2] to arrive at the 
effective nonlinear Schrodinger equation 



(4.2) idtU + -div, {{VlErn (Po)) ■ V,) ^ 



A(|z|' 



* \u\ )U, U\t=o = Uq. 



For /i < 0, existence of a smooth solution u S C{[0,Tc), H''), locally in time, can be 
proved along the same lines as in fTTl and hence, a result analogous to the one stated in 
Proposition |3.3| is straightforward. In the case /i > one can follow the arguments of |fT9]| , 
using a functional framework which is more intricate, however (the Sobolev spaces H'' are 
not sufficient but have to be intersected with weighted spaces), and we shall not do so 
here. In a similar spirit, stability in the sense of Theorem |3.4| follows from an adaptation 
of Lemma [331 which we leave out. 



4.2. Smooth kernels. If in ( 14. Il l the interaction kernel K(x) is a given smooth function, 
bounded as well as its derivative, then Ofc = 1 (corresponding/ormaZfy to the case /i = 0). 
Such a situation appears for example in [5], where 



Kix) 



(ai + a2|x|2 + a3|x|4)e-^ l-l 



-B2|j.|2 



with constants ai, 02, 03, 04 e M., A,B > 0. Resuming the above computations in this 
context, we find: 



K * IV'' 



e|2 



-d/2 



u t 



x-^- q{t) 



,P0 



K{C^e)\u{t,z-0\' 



i^(o) / |7.(i,^-c)rdc = i^(o)h(t)ii 
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due to mass conservation, along with an 0{^) convergence rate under suitable assump- 
tions. In particular, this shows that, as e ^ 0, the nonlinear effects become negligible. 
Indeed, in this case the envelope equation becomes 

idtu + ^div^ ((VfeSm (po)) • V^) u = K{Q)\\uo\\l2U, uit=o = "o- 
The right hand side involves a constant potential term, which can be gauged away via 

The remaining amplitude v{t,x) then solves a free Schrodinger equation with effective 
mass tensor V^i^m {Po)- 

Appendix A. Some useful algebraic identities 

For the derivation of the effective mass equation ( 11.81 ) we shall rely on several algebraic 
identities, which can be derived from Bloch's spectral problem (for more details see, e.g., 
lIU): First, taking the gradient w.r.t. to k of jl.H . we have 

(A.l) {H{k) - E„,) x™ + {H{k) - Era) Vfex™ - 

and, by taking the in L^(y)-scalar product with Xm, we obtain 

Since H{k) is self-adjoint, the last term is zero, thanks to (11. 2t . We infer 

(A.2) WkEjn{k) = (Xm, {-i^v + k) Xm)i2(y) ■ 

Differentiating ( IA.lt again, we have, for all j, ^ g {1, . . . ,d}: 

dl.k, {H{k) - E^) Xm + dk, {H{k) - E^) dk,_Xm + dk, {H{k) - d^^Xm 
+ {H{k)-E„,)dl,^Xm=0. 
Taking the scalar product with Xm, we have: 

dl^k,E,n{k) = 5ji + {{-idy^ + kj) dk.Xm + {-idy, + h) dk^Xm, Xm) ^2,y) 

(A.3) 

\^ki-^m^kj\m ^kj-^m^kt\m^\mf j^2^y) ' 

Appendix B. Adding an additional, slowly varying potential 

As a possible extension of our study, one might want to consider the case where the 
wave function is not only under the influence of the nonlinearity and the periodic potential, 
but also add an additional slowly varying external potential V{t, x), i.e. 

iedtr + = Vr (f ) + V{t,xW + e'^fiin^W- 

At least formally, this can be done by combining our analysis with the results given in ifToll : 
To this end, we define the semi-classical band Hamiltonian 

x) = E„,ik) + V{t, x), (fc, x)€Y* x 

and denote the corresponding semiclassical phase space trajectories by 

r m = Vfe£;„ , g(0) = go, 



(B.l) 



This system is the analogue of the classical Hamiltonian phase space flow, in the presence 
of an additional periodic potential Vpa-. In order to make sure that the system ( IB. II ) is 
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well-defined, it is sufficient to assume that Em{p{t)) is a simple eigenvalue {\Em{p{t)) ~ 
7^ for all n ^ m, t e M, fc G F*); see e.g. lITOl . where examples of such 
situations are given. 

The approximate solution under the form of a coherent state within the m-th Bloch band 
is then given by: 



with q[t) , obtained from dB.U . The highly oscillatory phase takes the form $,„ a-') = 
Smit) +p{t) ■ {x — q{t)), where Sm{t) G M is the (purely time-dependent) semi-classical 
action 



Note that <I>m simplifies to il.6\ in the case where V{t, x) = 0. In this case, the governing 
equation for the leading profile u{t, z) is found to be a nonlinear Schrodinger equation 
with time-dependent quadratic potential, time-dependent effective mass V'j,E,n{p{t)) and 
coupling constant \m{t), see |10| for more details. These features make it difficult to 
give sufficient conditions under which where the solution u{t, z) is global, i.e. Tc — oo. 
Indeed, the signature of WlE„i {p{t)) may change, and the existence of Strichartz estimates 
for the linear part is a non-trivial issue. Moreover, A„i (t) may also change sign, making 
the analysis even more delicate. 
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